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We present results of a combined density functional and many-body calculations for the electronic 
and magnetic properties of the defect-free digital ferromagnetic heterostructures obtained by doping 
GaAs with Cr and Mn. While local density approximation/(-|-U) predicts half-metallicity in these 
defect-free delta-doped heterostructures, we demonstrate that local many-body correlations cap- 
tured by Dynamical Mean Field Theory induce within the minority spin channel non-quasiparticle 
states just above Ef- As a consequence of the existence of these many-body states the half- metallic 
gap is closed and the carriers spin polarization is significantly reduced. Below the Fermi level the 
minority spin highest valence states are found to localize more on the GaAs layers being indepen- 
dent of the type of electronic correlations considered. Thus, our results confirm the confinement 
of carriers in these delta-doped heterostructures, having a spin-polarization that follow a different 
temperature dependence than magnetization. We suggest that polarized hot-electron photolumines- 
cence experiments might bring evidence for the existence of many-body states within the minority 
spin channel and their finite temperature behavior. 



INTRODUCTION 

Digital Magnetic Heterostructures (DMH) are semi- 
conductor heterostructures in which magnetic mono- 
layers are incorporated using the digital-alloy technique 
This technique was developed in order to com- 
bine the conventional semiconductors with magnetic ma- 
terials, with the desire to realize simultaneous band gap 
and magnetic engineering. Within the DMH technique, 
the quality of magnetic mono-layers is highly sensitive to 
growth conditions similarly to the case of diluted mag- 
netic semiconductors, were low temperature thermal an- 
nealing 0, 3] and codoping 0, @] has been extensively 
developed. The existing technological difficulties increase 
therefore the significance of first-principles investigation, 
which can predict many characteristics of DMH and in- 
dicate trends in their properties. 

For DMH materials, most first-principles investigations 
were done for GaAs/Mn 0, Q. The exchange interac- 
tions inside and between Mn mono-layers were analyzed, 
in particular the spacer dependent magnetic properties 
were recently discussed [§]. Independent of the spacer 
thickness the electronic structure was predicted to be 
half-metallic and using the technique of injection of free 
holes [10( intra and interlayer exchange coupling mech- 
anisms were discussed [§]. In addition, previous studies 
discussed transport properties, optimal concentration of 



Mn, possible effect of anti-site defects, etc., [7|, |8[ . A 
number of calculations were also performed for Mn doped 
Si, Ge, and Ga [ljjl3|. 



Properties of these heterostructures are expected to de- 
viate from the physical behavior presented by family of 
diluted magnetic semiconductors [14| . The most common 
position occupied by Mn within the (III,V) host GaAs 
diluted magnetic semiconductor is the Ga site. The sub- 
stitutional Mn contributes with its two 4s electrons in 
the crystal bonding similarly as the two Ga-4s electrons. 
Less favorable energetically is the interstitial Mn posi- 
tion and therefore is less common While in Mn 
doped GaAs semiconductors the system is described as 
"bulk", in the fabrication of GaAs/Mn-digital magnetic 
heterostructures, MnAs mono-layers were embedded into 
GaAs superlattices 15], so that the quasi-bidimensional 
geometry of the MnAs layer is preserved. Thus, the 
random substitution of positions (Ga occupied or inter- 
stitial) by the Mn ions in the case of diluted magnetic 
semiconductors shows up as an ordered replacement of a 
GaAs layers by MnAs layers in the case of DMH, with ev- 
ident consequences upon the critical temperature. DMH 
exhibits an unexpected Curie temperature dependence 
on concentration of Mn ions. At a fixed value of Mn 
concentration (fixed number of MnAs layers), the Curie 
temperature decrease with the thickness of the GaAs, 
but saturates at large thicknesses (= 50 GaAs mono- 
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layers) |15j . Electronic structure calculations were also 
reported in the literature for the case of i5-doping with 
Mn of GaAs/AlAs superlatices 0- In the defect free 
cases a half-metallic heterostructure was obtained. It was 
also demonstrated that in the presence of As-anti-sites, 
the half-metallic state was destroyed and a p-type metal- 
lic conductance was evidenced because of majority spin 
electrons. In addition the computed magnetic coupling 
between the Mn ions becomes larger in the presence of 
As anti-sites Q ■ To our best knowledge no results which 
include electronic correlations were reported for the case 
of Cr or Mn even in the defect-free cases. 

In this paper we present results of the first- 
principles electronic structure calculations within Lo- 
cal Spin Density Approximation (LSDA), LSDA+U and 
LSDA+Dynamical Mean Field Theory (DMFT) for the 
different defect free (MAs)i/(GaAs)7 digital heterostruc- 
tures, containing transition metal-arsenide mono-layers 
(M=Cr,Mn) introduced in the GaAs superstructure. The 
most important correlation effect in half-metals, the for- 
mation of many-body states within the gap can not be 
captured by the mean-field LSDA/+U-type calculations. 
Dynamical mean-field theory was demonstrated to be 
able to describe such effects [l6|, therefore is essential 
to compare directly the mean-field LSDA/+U and the 
DMFT results. We demonstrate that the half-metallic 
character is lost in <5-doped GaAs heterostrcutures con- 
taining Cr or Mn at finite temperatures and in the pres- 
ence of dynamic correlations. Correlation effects soften 
the magnetic properties of the ferromagnetic layers, in- 
duced states above Ep with tails that cross the Fermi 
level and reduce the carrier spin polarization. Within the 
minority spin channel below the Fermi level the character 
of states forming the band-edge remains not affected by 
electronic correlations and the states are mostly localized 
on GaAs layers confirming the experimentally observed 



confinement of carriers in these heterostructures 15, 17 



CRYSTAL STRUCTURE AND 
COMPUTATIONAL METHOD 

To compute the electronic structure for the digital 
magnetic heterostructures we used the standard rep- 
resentation of the tetragonal symmetry, space group 
P — Ami. The tetragonal super cell was obtained from 
the FCC unit cell by a rotation in the basal plane and 
an integer translation along the (001) direction. The 
new translation vectors are: a = a$/ \/2, b = ao/ V%, 
c = aoN, where ao is the bulk GaAs lattice parameter, 
and N = 4, and they correspond to a (GaAs)s super 
cell. Within this super cell a GaAs mono-layer is re- 
placed by MAs, such that the obtained DMH has the 
unit cell formula (MAs)i/(GaAs)2N-i- For GaAs which 
has an open structure, a close packed structure is ob- 
tained by including empty-spheres. Such empty potential 



wells were used also in the present heterostructure geom- 
etry. We used identical muffin-tin spheres for all atoms 
in the unit cell having the average Wigner-Seitz radius 
2. 62852a. u. corresponding to the experimental GaAs lat- 
tice parameter (ao = 5. 65 A) . The basis used for the 
self-consistent calculations contains the spd-partial waves 
for transition metals, sp(d)-partial waves for Ga, As and 
s(p)-partial waves for the empty spheres E. (I) means 
that the Z-partial waves are downfolded within the self 
consistent calculations. Selfconsistency was performed 
for 105 irreducible fc-points and total energy convergence 
was achieved with an accuracy of 10~ 6 Ry. Results for 
the density of states, magnetic moments and spin polar- 
izations are presented in the next section. 

Correlation effects in the valence 3d orbitals are in- 
cluded via an on-site electron-electron interaction in the 

form J Ei{m,(r} ^mm'm"m'" c j rtw c i m v' c im"VC»ro"<r. The 



interaction is treated in the framework of DMFT [18M20j . 
with a spin-polarized T-matrix Fluctuation Exchange 
(SPTF) type of impurity solver 21] implemented within 
the EMTO basis set [ll HI]. The Coulomb matrix ele- 
ments Umm'm"m"i are expressed in the usual way [24| in 
terms of three Kanamori parameters U, U' = U — 2 J 
and J. Here, Ci mCT /c| mcr destroys/creates an electron 
with spin a on orbital m on lattice site i. The SPTF 
approximation is a multiband spin-polarized generaliza- 
tion of the fluctuation exchange approximation (FLEX) 
[25l [26} , but with a different treatment of particle- hole 
(PH) and particle-particle (PP) channels. The particle- 
particle (PP) channel is described by a T-matrix ap- 
proach 27J, l28j giving a renormalization of the effective 
interaction. This effective interaction is used explicitly 
in the particle-hole channel. Justifications, further de- 
velopments and details of this scheme can be found in 



Ref. 21 



For the case of half-metallic ferromagncts it was 
demonstrated [l6| by model as well as realistic electronic 
structure calculations that many-body effects are cru- 
cial for half-metals: they produce states with tails that 
cross the Fermi level so that the gap is closed and half- 
metallicity is lost (29H33j . The origin of these many- 
body non-quasiparticle (NQP) states is connected with 
"spin-polaron" processes: the spin-down low-energy elec- 
tron excitations, which are forbidden for the HMF in 
the one-particle picture, turn out to be possible as su- 
perpositions of spin-up electron excitations and virtual 
magnons [3, H3] • Spin-polaron processes are described 
within the SPTF approach by the fluctuation potential 
matrix W aa (iuj) with a — ±, defined in a similar way as 
in the spin-polarized FLEX approximation 261 ]: 



W ++ {iuj) W+-{iuj) 
W-+\iJ) W (iu) 



(1) 



The essential feature here is that potential (JTJ) is a 
complex energy dependent matrix in spin pace with off- 



3 



diagonal elements: 



W° 



(2) 



where U m represents the bare vertex matrix correspond- 
ing to the transverse magnetic channel, x <T ' _!T (iw) is an 
effective transverse susceptibility matrix and Xo :-CT (* a -0 IS 
the bare transverse susceptibility j26| . iuj are fermionic 
Matsubara frequencies and (m) corresponds to the mag- 



netic interaction channel 



The local Green func- 



tions as well as the electronic self-energies are spin diago- 
nal for collinear magnetic configurations. In this approx- 
imation the electronic self-energy is calculated in terms 
of the effective interactions in various channels. The 
particle-particle contribution to the self-energy was com- 
bined with the Hartree-Fock and the second order contri- 
butions [i^l . To ensure a physical transparent description 
the combined particle-particle self-energy is presented by 
its a Hartree TS TH \ioS) and Fock ^ TF \iuj) types of con- 
tributions: E(iw) = ^ TH \iu) + ^ TF \iuj) + ^P h \iuj) 
where the particle-hole contribution E^) reads: 



J l2o 



34cr' 



(3) 



Since the static contribution from correlations is already 
included in the local spin-density approximation (LSDA), 
so-called "double counted" terms must be subtracted. In 
other words, those parts of the DFT expression for the to- 
tal energy that correspond to the interaction included in 
the Hubbard Hamiltonian has to be subtracted. Several 
double-counting schemes has been proposed and used for 
specific materials applications 3^, 3(| . It has been recog- 
nized that the results of the LSDA+U calculations may 
depend crucially on the choice of the scheme used [371 ]. 
In particular for the class of moderately correlated met- 
als as one may consider the half-metallic ferromagnets 
in the LSDA+U calculations performed here we used 
the "around mean-field" [37I ] double-counting scheme. 
Similarly to the LSDA+U in the DMFT approach the 
Hubbard Hamiltonian represents the underlying physics, 
while in contrary to the mean-field LSDA+U solution the 
time dependent dynamics are considered, while spatial 
fluctuations are neglected. Obviously double counting 
correction has to be considered also. To achieve this, we 
replace Y> a {E) with T, a (E) — E ff (0) [II] in all equations 
of the DMFT procedure [19( . Physically, this is related 
to the fact that DMFT only adds dynamical correlations 
to the LSDA result. For this reason, it is believed that 
this kind of double-counting subtraction "£(0)" is more 
appropriate for a DMFT treatment of metals than the 
alternative static "Hartree-Fock" (HF) subtraction [37| . 



RESULTS AND DISCUSSION 

When the mono-layer of MAs, with M=Cr, Mn is in- 
troduced in the GaAs super cell, the symmetry lower- 



ing from cubic to tetragonal take place and the 3c? — ti g 
states split further into a double degenerate d xz and d yz 
situated at lower energies and the non-degenerate d xy . 
At higher energies the d^ z 2_i and d x 2_ y 2 are situated 
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FIG. 1: (color online) Spin-resolved total DOS of 
(MAs) 1/ (GaAs) 7, with M=Cr and Mn. The half-metallic 
character is evidenced by the simultaneous presence of a 
metallic majority-spin and an insulating (semiconducting) 
minority-spin DOS obtained within LSDA (dashed black-line) 
and LSDA+U (dot-dashed blue line) calculations. Dynami- 
cal correlations captured by DMFT (solid red-line) introduce 
many-body states within the gap that destroy half-metallicity. 



24 |3Jj. Magnetic exchange interaction split the spin 
channels and contribute additionally to the crystal field 
in opening the minority spin gap. In Figfl] we present 
spin resolved, total density of states (DOS), for the DMH 
structures (MAs) 1/ (GaAs) 7, M=Cr, Mn. In both mi- 
nority and majority spin channels at energy of about 
— 12eV (not shown in FigsQ]) the As-4s and Ga-4s are 
present. In the energy range —8 to — 2eU As— p and 
Ga— p states are dominant characters in the p — d bond- 
ing orbitals, while in the vicinity of the Fermi level in 
the energy range of —2 to 2eU the dominant 3d transi- 
tion metal states are present (see also orbtial projected 
DOS of FigH]). At the LSDA level the occupied minority 
spin channel for both compounds is similar, because in 
this energy range the transition metal d-states are weakly 
present within the p — d bonding orbital dominated by As 
and Ga-p states. The Fermi level is situates for minor- 
ity spin electrons within the gap, closer to the bottom of 
the conduction band. The values of the energy gaps are 
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&LSDA ^ 51eV and ^lsda ^ 0A2 eV. The difference 
between the Fermi level and the bottom of the conduction 
band is somewhat larger for the DMH with Cr, in com- 
parison with the value obtained in the Mn-based DMH. 
Given the fact that the top of minority-spin valence band 
is situated at similar energies in these heterostructures, 
the magnitude of the gap is therefore determined by the 
energy distance between the Fermi level and the bottom 
of conduction band. 

A more detailed analysis of the density of states can 
be made on the base of the orbital projected densities 
of states Fig. H Here the LSDA and the LSDA+U re- 
sults for the correlated transition metal atoms are com- 
pared. LSDA+U calculations were performed for differ- 
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FIG. 2: (color online) Orbital projected spin polarized DOS 
for (MAs)i/(GaAs) 7 , with M=Cr (left column) and Mn(right 
column). The LSD A/LSD A+U results are shown with black- 
dashed/blue-solid lines. 

ent values of the average Coulomb interaction U = 2, 3, 4 
and 5eU taking the same value for the exchange param- 
eter J = 0.9eV. In Fig. [2] we show the results for 
U=3eV and J=0.9eV, typical values for transition met- 
als 24, [4(|. The left column represents the results for 



the densities for the compound containing Mn. In each 
pannel the upper/lower part contains the results for the 
majority /minority spins. Although the similarities of the 
orbital projected DOS are visible, we note the general 
tendency of the LSDA+U namely to enlarge the minority 
spin gap A^ s r DA+u = 0.85eU and A^ A+U 0.81eU. 
An interesting feature is that most of the spectral weight 
at the Fermi level is provided by the 3d x 2_ y 2 orbital for 
both cases of Cr and Mn, however the larger value is 
obtained with the Cr substitution. An almost negligible 
spectral weight is obtained for Sd^-orbitals. As expected 
from this approach transition metal d-states are further 
pushed apart from the Fermi level. In particular for the 
Cr heterostructure within the majority spin channel the 
LSDA peaks (see Figs. II 1211 located at about -leV are 
shifted to lower energies at -1.5eV. The occupied ma- 
jority spin Mn states are within the bonding band at 
-3eV are also shifted to lower energies at -4eV. Within 
the minority spin channel for all compounds the group 
of states situated around leV above Ep are shifted ac- 
cordingly to higher energies (see Fig. [TJ [2]). For the 
other values of U and J a similar tendency of DOS was 
obtained for both heterostructures. In the followings we 
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FIG. 3: (color online) LSDA (black-dashed) and LSDA+U 
(blue-solid) As (left column), Ga(right column) spin polarized 
DOS for (CrAs)i/(GaAs)7. The positions in the unit cell 
are specified within the figures. The U=3eV and J=0.9eV 
parameters were used in the LSDA+U calculation. 



the compound with Cr, while the right column shows plot the total As and Ga spin resolved density of states for 
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(MAs)i/(GaAs) 7 , M=Cr,Mn. Fig|3] presents the results 
for Ci-5 doping. The left/right columns show the As/Ga- 
DOS for the different layers. As one can see in the major- 
ity spin channel the most important contribution around 
the Fermi level comes from the arsenic and gallium atoms 
situated at As : a/2(l, 0, c/8a), a/2(0, 1, — c/8a) and Ga : 
a/2(l, 1, c/8a). These are in fact the positions in the im- 
mediate vicinity of the transition metal Cr As-layer. The 
similar trend is visible for the Mn-5 doped heterostruc- 
ture in FigJU In addition we note that as expected at the 
Fermi level the dominant states are As-p states. The gen- 
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FIG. 4: (color online) LSDA (black-dashed) and LSDA+U 
(blue-solid) As/Ga spin polarized DOS for (MnAs)i/(GaAs) 7 
are shown on the left /right columns. The positions in the unit 
cell are specified within the figures. The U=3eV and J=0.9eV 
parameters were used in the LSDA+U calculation. 

eral picture resulting from the extended analysis (by in- 
specting the atom and orbital resolved densities of states 
-FigsH |3] and |H ) of the LSDA and LSDA+U results is 
that the £-doped heterostructure presents a half-metallic 
behavior with a minority spin gap. 

LSDA+DMFT calculations were performed for the 
same values of U = 3eV and J — 0.9eV and temper- 
atures T = 30, 40, 50K . In Fig. we present the DMFT 
results for the atom resolved total density of states in a 
small energy window around Ep. Further we show the 
3d-orbitals contribution into the main peaks of Cr/Mn 
transition metals. The insets of Fig. [5] shows the results 



of non-magnetic calculations, from where we can roughly 
estimate the bandwidth which is of about 6eV. Based on 
the analysis of the values of the bandwidth, orbital oc- 
cupations and the magnitude of U, one expects for the 
heterostructures in discussion, a behavior typical for a 
medium correlated metal. 
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FIG. 5: (color online) LSDA+DMFT atom resolved spin po- 
larized DOS for (MAs)i/(GaAs) 7 , with M=Cr (upper panel) 
and Mn(lower panel). The results presented were obtained for 
the following parameters: U M n/Cr = 3eV, JmhJCt — 0.9eF 
and T = 50AT. For the correlated atoms the main orbital 
contributions are also shown (solid red line). Average total 
As (dashed blue) and average total Ga (green dote dashed) 
contributions are seen also. Above Ef, many-body induced 
states with predominant d xy character determines the closure 
of the half-metallic gap. Insets show the orbital resolved non- 
spin polarized Cr/Mn-LDA density of states results. 



As one can see the dynamic correlations captured by 
DMFT gave a completely different picture for the or- 
bital distributions around the Fermi level. In the case of 
Cr significant majority spin density of states is obtained 
for the Sd^y-orbitals in contradiction with the LSDA/+U 
results were these orbitals are shown to give a negligible 
contribution see Fig. 03 Above Ep + O.SeV, in the unoc- 
cupied part the majority spin 3d electrons have a signif- 
icant 3d x 2_ y 2 and at higher energies (not shown) 3d z 2_ 1 
character. In the minority spin channel from Ep — 0.5eV 
no significant spectral weight is present, while above Ep 
the many body induced states have a predominant 2>d xy 
character. At higher energies a predominant 3c? z 2_ 1 char- 
acter is present. 
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In the case of Mn-S doped structure, in the major- 
ity spin channel As-p is dominant at the Fermi level, 
the Mn-3d x a_„2 contribution is seen around 0.3eV below 
the Fermi level. It is interesting to note that the posi- 
tion of this orbitals are not changed with respect to the 
LSDA picture. Contrary to the LSDA minority spin Mn- 
3d x 2_ y 2 have a larger spectral weight below Ep, while 
above Ep 3(4 y -states are seen to cross the Fermi level. 
The general picture obtained by including dynamic cor- 
relations is that the gap within the minority spin channel 
is closed and the half-metallicity disappears. For both Cr 
and Mn compounds closure of the gap is due to the pres- 
ence of many-body states with 3d xy character only seen 
in DMFT. 

The analysis of the orbital resolved density of states 
can be correlated with the behavior of the self-energy 
seen in Fig|6l where the atom and orbital resolved imag- 
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FIG. 6: (color online) Orbital resolved imaginary part of mi- 
nority/majority spin self-energy for Cr (red solid/dashed) and 
Mn (green solid/dashed) lines computed for UMn/o=3eV 
and J=0.9eV and T=50K. The hump at around 0.2eV vis- 
ible for d xy orbital signals the departure from the expected 
Fermi liquid behavior. 

inary part of self-energy is shown. The self-energies for 



d za: -orbital were excluded since they are identical due to 
the cell symmetry with the d yz ones. We observe that the 
imaginary part of the self-energy has a rather symmetric 
energy dependence around the Fermi level, with a normal 
Fermi-liquid-type behavior ImY^ (E)c r /Mn ^ {E—Ep) 2 . 
The minority spin ImYr r {E)cr / Mn shows a significant 
increase just above the Fermi level which is more pro- 
nounced for the d xv orbitals. 




E-E F (eV) 

FIG. 7: (color online) Orbital resolved real part of minor- 
ity/majority spin self-energy for Cr (red solid/dashed) and 
Mn (green solid/dashed) lines computed for UMn/cv=3eV 
and J=0.9eV and T=50K. A solid black line with a unitary 
negative slope is also shown. 

In Fig. [7] we plot the real part of the orbital re- 
solved majority and minority self-energies. The left/right 
columns represent the results for Cr/Mn doped com- 
pounds. The figure also displays a line with a unitary 
negative slope as a guide for the analysis of the energy 
dependence of the real part of the self-energies around 
Ep. As one can see for both spin directions, the self- 
energy has a negative slope <9£(w)/<9w at the Fermi en- 
ergy, which confirms that the quasiparticle weight Z = 
(1 — dH{iS) I duj)^ 1 is reduced by correlations. For both 
Cr and Mn for majority spins <9£(u;) / du is clearly less 



7 



than unity. Contrary for the minority Cr/Mn-dxy spins 
dT, xy (uj)/doj is larger than unity above the Fermi level 
(within our approximation, we cannot determine with 
sufficient accuracy), suggesting the nonquasiparticle na- 
ture of the minority-spin states within the gap. We have 
checked that the above behavior of the real and imagi- 
nary parts of the self-energy does not change significantly 
in the range of investigated temperatures. For smaller U 
values the slope is reduced slightly although, the same 
peculiar behavior is seen for the real and imaginary part 
of d^-orbital self-energy. Therefore we can conclude that 
the many-body induced NQP states within the minority 
spin gap situated just above the Fermi level, are mainly 
determined by the transition metal d xy orbitals. As the 
closure of the minority spin gap happens by the extension 
of NQP state tails from above, towards the Fermi level 
and in order to have a complete picture it is important 
to understand the effect of correlations upon the minor- 
ity spin states situated below Ep. This can be studied 
by looking at the band edges within the minority spin 
channel. In particular we investigate the character of the 
electronic states below and above Ep with or without in- 
cluding correlations. As discussed previously, above the 
the Fermi level transition metal 3d xy states provide the 
main character that developed many-body NQP states 
who's tail cross Ep. Below Ep the dominant states have 
As-p character. Therefore, we looked at contributions of 
As states originating from different As layers to the top of 
minority spin valence band. The As layers are grouped 
into four classes, each having the distance za s — f ff , 
where j = ±1, ±3, ±5, ±7 with respect to the transition 
metal atom located in origin. The results of this analy- 
sis is presented in table |TJ As the distance between the 





(CrAs)i/(GaAs) 7 


(MnAs)i/(GaAs) 7 


Px Py 

% % 


Px Py 
% % 


Asl:f(0,l,£) 

Asl:f(l,0,-£) 

As2:f(0,l,f§) 

As2:f(l,0,-ff) 

As3:|(0,l,|f) 

As3:|(l,0,-|f) 

As4:|(0,l,|f) 

As4:f(l,0,-£) 


10.14 16.66 
27.57 22.22 
30.43 30.55 
32.88 30.55 


10.11 17.23 
27.03 20.43 
30.18 31.70 
32.66 30.64 



TABLE I: Relative contribution (in percent) of the As-p states 
to the valence band maxima from each As-layer obtained 
within the LSDA calculation. Local correlation within the 
mean-field LSDA+U or DMFT does not change significantly 
these values. 

transition metal mono-layer and As mono-layer increases, 
the As-(p x , p y ) character continue to increase. In the 
same time, the difference between values from different 
mono-layers decreases. These results suggest that the As 
mono-layer situated closer to the transition metal has its 
electrons confined by the p — d hybridization. At larger 



distances, hybridization decreases, and more p-character 
is available to form the top of the valence band. One can 
note that p x and p y contributions arc different as a con- 
sequence of symmetry lowering present in the tetragonal 
geometry. These results are in agreement with previous 
calculations obtained using the macroscopic average of 
the total DFT potential 01111 Ell- 
in the followings we comment the comparison of the 
results obtained within the three different approaches to 
the electronic structure calculations within the DFT. A 
well known problem of the LSDA and its relative GGA is 
the underestimation of the band gaps (for many semicon- 
ductors) in addition to the missing electronic correlation 
effects. For the heterostructures in discussion we present 
results obtained with different methods that consider cor- 
relations effects: a simplified mean field LSDA+U and a 
more complex LSDA+DMFT approach. The LSDA+U 
scheme seeks to correct the problems of LSDA by adding 
a repulsive potential +U, obtained from a mean field 
decoupling of the local Coulomb interaction. Its con- 
sequence is that the occupied/unoccupied 3d levels are 
shifted to lower/higher energies as seen in Figs. I2I3I4I 
The ground state picture resulting from both this meth- 
ods is a half-metallic ferromagnetic state. 

On contrary to the mean-field methods, DMFT cap- 
tures correctly the dynamics at low as well as high en- 
ergies. This approach describes properly the quasipar- 
ticle formation around Fermi level simultaneously with 
the lower and upper Hubbard bands at higher energies 
[Lsl Il9j |. Dynamical correlations captured by DMFT 
change completely the physical picture of these com- 
pounds. For the majority spin states a slight spectral 
weight redistribution takes place, however this does not 
affect the value at the Fermi level which is similar to 
that obtained within LSDA (see Fig. []}. For minority 
spins, although the position of the top of the minority 
spin valence band is slightly changed including correla- 
tions, the As-p contributions from different As layers are 
similar in values to those obtained within the LSDA cal- 
culations, see Tab. |U This demonstrates that the As- 
p character distribution within the top of the minority 
spin valence band is not sensitive to correlation effects. 
The conclusion therefore is that the many-body inter- 
actions affects the half-metallic band structure predom- 
inantly within the minority spin channel, above Ep, by 
introducing NQP states. The presence of these many- 
body states reduce significantly the ideal 100% carrier 
spin polarization predicted by LSDA(+U). In order to 
discuss quantitatively the conduction electrons spin po- 
larization we use the common definition [43TJ45l |: P n = 
(N t {E F )v^(E F ) - Ni (Ep)v™ (Ep ) ) / (N^(Ep)v™ (Ep ) + 
N±(Ep)v™(Ep)), where N and v are the density of states 
and the velocity at Ep. The values of n determine what 
type of experimental measurements that can access the 
polarization. The Pq corresponds to the spin-resolved 
photoemission experiments, while the higher orders cor- 
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FIG. 8: (color online) LSDA+DMFT values for the fi- 
nite temperature magnetization and spin polarization of 
(MAs)i/(GaAs) 7 , M=Cr, Mn compounds. The LSDA val- 
ues are taken as T=0K results. 



responds to spin polarizations measured with point con- 
tact Andreev reflection in the ballistic (n=l) or diffu- 
sive regimes (n=2) [HI, 0|. In Fig. |8j we show the re- 
sults for the polarization obtained using the simplified 
Pdmft{Ef ,T) = P n —o, where the velocities for spin up 
and down are droped. One can see a clear distinction 
between the different temperature behavior of magneti- 
zation and polarization, similar to many potential half- 
metallic ferromagnets [H EI 0, E3] . 

In table [TT| we present the values of the finite tempera- 
ture magnetizations, the values of polarization being eas- 
ily readable from the lower panel of Fig. [8j For the 





Mlsda 


Mlsda+u 

fJ-B 


Mdmft 


30K 40K 50K 


(CrAs)i/(GaAs) 7 
(MnAs)i/(GaAs) 7 


3.00 
4.00 


3.00 
4.00 


2.97 2.82 2.76 
3.80 3.67 3.52 



TABLE II: Integer magnetic moments values were obtained 
within LSDA(+U). The finite temperature DMFT magnetic 
moments demonstrate departure from the half-metallic val- 
ues. 

Mn-S doped GaAs the temperature dependence of rem- 
nant magnetization of samples that consist of 100 period 
of GaAs(10 mono-layers)/MnAs(0.5 mono-layer) showed 
a quasi-linear dependence and a critical temperature of 
about T c « 50K [l5(. More recent measurements on the 
Mn-<5-doped DFH include polarized neutron reflectome- 
try and magnetometry performed on samples having 20, 
50 and 100 mono-layer GaAs spacers in the tempera- 



ture range from 30 to 40K [48[ . This analysis allowed to 
identify two types of contribution to the magnetization 
of Mn-doped DMH: one type that has the superlatticc 
periodicity and the second one resulting from random 
clustering effects suggested by the discontinuous nature 
of the quasi- two-dimensional MnAs 48] . A direct com- 



parison with these experimental results are not possible 
for the present theoretical calculations, due to the fact 
that the supercells considered here contain a considerably 
smaller number of layers. Nevertheless certain qualitative 
agreements are obvious. As one can see in Fig. [5] the 
present LSDA+DMFT calculation captures a qualita- 
tive quasi-linear temperature dependence of the reduced 
M(T) /Mlsda magnetization also for the smaller num- 
bers of mono- layers in discussion, a behaviour that is seen 
in the experiment. In addition the quasi-two-dimensional 
nature of the problem is reflected in the present calcu- 
lations by the separation of electrons and holes among 
different layers (see Tab. Hand the discussion below). A 
major difference is that for the current supercell geome- 
try our estimate for the Curie temperature is around four 
times larger. 

Let us further comment upon the recent polarization 
measurements of spin carriers in Mn-<5 doped GaAs com- 
pounds using hot electron photoluminescence (HPL) [13] ■ 
The hot-electron photoluminescent spectra were inter- 
preted based on the previous analysis performed in di- 
luted magnetic semiconductors 49]. According to this 



analysis in Mn doped GaAs in a wide range of Mn content 
the valence band holes predominantly occupy the Mn ac- 
ceptor impurity band. It was suggested that the polar- 
ization of the valence band holes is proportional to the 
sample magnetization 



17 



49(. To interpret the data on 
digital Mn/GaAs heterostructures, a combined contribu- 
tion from the (5-doped and the unintentionally doped in- 
terlayer region was considered as a superposition of both 
contributions. It was shown that carriers experience a 
strong exchange interaction with the ferromagnetic in- 
terfaces, which rapidly decease with increasing distance 
with respect to the ferromagnetic Mn-<5 layer [17], 

As one can see the electronic correlations play an im- 
portant role in the properties of the defect-free DMH 
presented above. In the presence of dynamical many- 
body correlations the top of the valence band does not 
change significantly, however, NQP states appear in the 
minority spin channel, FiglT]and[5l just above the Fermi 
level. At finite temperatures the tails of the NQP state 
cross the Fermi level and the half-metallicity is lost. In 
addition we have shown that carriers spin polarization 
and magnetization follow a different temperature behav- 
ior. We expect that optically excited radiative recom- 
bination of electrons and holes within the minority spin 
channel would help to demonstrate the presence of the 
many-body effects in particular the existence of non- 
quasiparticle states in these materials. Hot-electron pho- 
toluminescence spectroscopy was recently used to dis- 



9 



cuss hole spin polarization in diluted magnetic semicon- 
ductors [5y]. The HPL circular polarization under cir- 
cularly polarized excitation provides detailed informa- 
tion on spin-relaxation mechanism [49l . 51|. and allows 
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a finite temperature characterization of the HPL polar- 
ization [H(J. We believe that such experiments would 
demonstrate a different finite temperature behavior of 
magnetization and polarization as qualitatively demon- 
strated in Fig. |U 



CONCLUSIONS 

The technology based on digital alloy technique is ex- 
pected to produce heterostructures with an enhanced Tc 
because of locally confined high transition metal concen- 
tration, and a better control of the ferromagnetic prop- 
erties by controlling the GaAs layer thickness. The ex- 
pected enhancement of the Curie temperature has not yet 
been realized, instead, it was found that magnetic prop- 
erties depend strongly on inter layer thicknesses [la, |48| . 
Our results of the electronic structure calculations based 
on LSDA predict that the defect free delta doped dig- 
ital magnetic heterostructures are half-metals in agree- 
ment with previous calculations Q- The same conclu- 
sion is given by computations that are using a mean field 
LSDA+U approach. On contrary including many-body 
correlations captured by DMFT the half-metallicity is 
lost. The computed finite temperature magnetization fol- 
lows an almost linear temperature dependence similarly 
to the experimental measurements performed however for 
larger GaAs(10 mono-layers) /MnAs(0.5 mono-layer) het- 



erostructures [15|, |48j . In comparing the critical temper- 



atures, our results overestimate by a factor of four the 
measured values. Certainly additional work is necessary 
in this direction in particular the extension of the the- 
oretical approach to include dynamical correlations be- 
yond locality. A particular importance in applications of 
DMH in spintronic devices is the effect of the ferromag- 
netic layers on the spin polarization of the carriers. Our 
results demonstrate that magnetization and polarization 
follows a different temperature dependence. We suggest 
that such effect might be captured by hot-electron pho- 
toluminescence spectroscopy. The quantitative analysis 
of the intensity of recombination radiation for transition 
between the electrons and holes within the minority spin 
channel in the presence of electronic correlations is in 
progress. 
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